SHilHHF »s-n«a .--.arA ns r 

sSHSr* r;r ~_ :r - 

corresponds to radially constant displacement velocity Th^ av r*i * , 

components of the vortex sheet motloS are then given ^ S "' rl 

“axial = v ' cos2 ♦. ( 


V* cos 4> s sin 0 


MOTION OF THE ENTIRE SLIPSTREAM 

fh.i/h^! 1 ' ix an a PP endi * to reference I, pointed out that the slipstream 
fluid between the vortex sheets moves at a fraction F of the sheet velocit? 

si; 


4 . b /r + i 
2 3J 




Here A Is the advance ratio 


A = V/fiR = (V/nD)/7r 




w . . = F 
axial 


W swirl = F v ’ cos K s '" t 


F = (2/tt) cos -1 (e -f ) 
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THE RADIAL CIRCULATION DISTRIBUTION 


slipstream motlon'l^fjJnd^ 0 ^ this minimum induced loss 

equal to the total vort Id tv trailed hv ! h ? T ? I °1 ab ° U+ a sl 'Ps+ream tube 
and introducing -light loading- app^x imIJtoni? ^ COrrespond '"3 radius. 


BT - 2ir r s F v’ cos <j> sin <J> 


r s = r 


*S = ♦ 


( 6 ) 

(7) 

( 8 ) 


(9) 


<t> = tan *(V/nr) 

The resulting circulation function is conveniently. written in a normalized form 

2 


_ F x‘ 


x 2 + I 


= G 


( 10 ) 


(G for Goldstein or Glauert) where 


x = Qr/V = Cr/R)/X 


(II) 


Goldstein, in his doctor's 


Equation (o seems too simple to be true 

operating 6 ^ & however, for prope! iers 

are nearly f,a?" paraHeT anS c?oJ! man ' , J b lades ' “here the vortex sheets 

circulation functions are'compared in figure 3 It shou I d ’ : Be+z . a J d 60 Ostein 
radial plot of G is identir^l wi+h .? T® , should be noted that a 

axial slipstream velocity ( increment ) r ?o l ?he P d?I P ?L + m:nt ra v + i?ocltr e 


DETERMINATION OF THE DISPLACEMENT VELOCITY 


disc ^oadTngHft^ust ^or n torque ) ^hy^eso ! + | S d 1 sp 1 a ^ emen+ velocity to the 
nal components: ^ v,n 9 oukowsky’s law into two orthogo- 


(jrjl m 00r (I - a') TB (12a) 

(r dr^L PV (l + a) rB <l2b) 

Here the sebscript " L » „eans that only the ,i,t forces are being considered. 
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I I ft^ng a ^fnes? re Retain'f , ng l the C *nght * load^n^as 2 °V h6 ,nduCed velocit V a+ the 

velocities to be half the vortex shee?2e?oc l tie^ln^ S° d I 3 *’" 9 | he induced 
we obtain: cities in the developed slipstream. 


= ' - + s s'" » s 5 y fe, -j 


X + I 


, - 1 V 2 ^ | f v ») 

3 - o- -w- cos d> = J- _ 

2 V v s 2 |V I 2 


x + I 

The circulation is given by the Betz-Prandtl approximation: 

r - [*&)* 

The blade profile drag contributions to thrust and torque are given by 

fell ■ 

dr D 


(13a) 

(13b) 


(10 restated) 


f°] dL . „ 

■ (rj d? s,n * = - 

fi dQ| _ . fD) dL . o. 

[r drj D " + [lJ d? COS * = + 


f— ) 1 

(drj L x 

(14a) 

Ml x 

[r drj L * 

(14b) 


The radial gradients of thrust coefficient T = ?T/r,v2_ D 2 j 

dent. P c = 2P/pv3 7rR 2, fInal|y may ^ wNtlehs ^ ^ ^ P ° WSr COeffl - 


dT c dI | 


dl 


2 _2 


d? d £ ^ “ d£ C 


d p c dJ 


dC dT dT 


dd 2 2 
+ -rJ- C 


(15a) 


(15b) 


where S = r/R, ? = v'/V, and 


dI . 

dl 


= 4 5 G 


dT 2 

d? 


= 256 1 


dJ 


I 1 - 




(16) 

(17) 

(18) 
(19) 
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four + !n+egrals l |’ an<J 19 ^ ai ? be k ? u [" erIcal ly Integrated radially to give 
sumed radial di sir? but ion 'of "prof i 7 e ' D/L^at io T* B ^ +h ® Pr ®“ 

-tic is then easily found wifh these filing: 


C = 


21 . 


j; 


4 T 


c 2 


(thrust) 


(20a) 


? = 



(power) 


(20b) 


attack of an el I ipticaM y loadid\ing^C | l /u(b 2 /Sl® rpar+ ° f +he induced an 9 'e of 


power coefficient, p£, 'and^he^i'sptacement 9 ^^'^-^ P ° W ? r ' one calculates the 
20b; the thrust coefficient and the effir- velocl+ y ra+ ‘o. C. from equation 
and n = T./P„, resDecHve?u ?hI h !,^l!^! nCy are +he " by T c ■ r r , .2 


and o ■ T " yp c; o b®apect i ve I y . The a I ternate "procedure^whe^the^hrust 1 ! s~sj5ec- 


i f ied, is 


or,, + ^° r m ° dera+e| y loaded propellers operating 
equatjons 20a or 20b may give values of C which 

IJ! S . ca f. a second approximation of the radi< 
coefficient is given by 


at low advance ratios, 
are large compared to A. 

I gradients of thrust and 


power 


where 


and 


dT , . 

"d r “ 4 C A G ( t ) ( cos 4 > - £ s i n d>j 

dP f w 

HT = 4 C ^ G [yj ( s ' n <}> + ° cos J 


<j) = tan 



+ 


\ 

C 

2 , 


(21 ) 
( 22 ) 


( 23 ) 


W 

V 



y C cos <p 


( 24 ) 


of 


a E nd Ua ^°appropr?etftr?hryal« 7?"* f ° better Values 

Fol lowing Theodorsen fret. 4 , . „ ^p^lmation 

r-wr-e — _ X • II \ ... 


n which G(A,B) is recalculated with 
to account for slipstream distortion 


it . . “ ■'r.se, u* i mar ion 

vortex advance ratio", A y = Ad+5/2), 


determination of- the propeller geometry 


coef f lctent°f or ' the^eq u I red ^c i rcu I at i on : COndTo * * e< * b V choice of lift 

1 , w 2 

2 P y cc £ =p wr=p W 


2t rVv 1 


This can be written as - = 4lT * 

R 


B (W/V) c. “ B 


&n 

G 


C_ ~ 4ttX 

B AT 

/x t 


(G) 

(25) 

? 

C * 

(26) 


th[ckness-to-chord m ratlos alTnV'LV 1° s ! rUCtural “"faints on 

also they must be consistent vith the D/L°rat os’thnt'!! "'"If 5 ot outer r ad i i; 
Some consideration must be aiven to f t,os +hat have been used to find £. 

a propeller designed ?or cru se cln bl^THT^' 0 " 5 35 WeM; for a * a *P>e, 
cient increases at inner rad i than It f pec+ed T deVelo P lar W ,ift coeffi- 
advance ratios, as in climbing flight! U ^ rad " Whe " ' + ' S opera+ed a + lower 

such as the 'c'l'ark^^^Constder i'n^th^'l'arqe ^th^' b ° +t0m aJrfoil sec +i°ru 
structurally at the inner radii 9 a nd th^ ? +hlc J ness ~ to - ch ord ratios needed 

:ener a To e ^? a r ?r» r !;b n rLIe, TO ' cln telgn^h^ro 
the propeller wl I , ^ ti!,. cese P 

P 


r , 

geometric . [ | ) 

Diameter I + 2 C I » ot = 0 


(27: 


Modern computational airfoil theorv (ref 5 } A _ 

for dark Y airfoils of varying thicknesi-to-cfTd X 7s g i ll™' ^ 


C l - 0.062 + 4.21 (t/c) + 0.0971 ot° ; 0.07 < 1 < o. 19 
when they are operated at a Reynolds number of I * and a Ma ch number of 0. 

through%he“i:openerdis C S at f vanishlnllv"? d n'* 0 ™ fl °" at fli 9ht velocity V 
realistic assumption for propeMers 1? -'“'f; ? f Thls ' = "°+ a 
are often quite large compared to the rec ' b . dri ve P lst on engines whicl 

distribution, averaged aroSnd Ihl p^pe?^? 6 aT^dl ' f ^ a *' a ' ve, ° ci ^ 

<t is customary to "depitch the propeller" (£(E)<|} sf that V TT®" by UV ’ 
given by K p 1 * so that the blade angle Is 


8 = tan - * 


f (■ + §) 


+ a 


(28 


The S performance f consequences r of n Dron e | pres ^ r ' b ® d circulation function, 
in the next section. P pel ler_fusela 9e interaction are considered 


290 




ORIGINAL PAGE 1$ 
OF POOR QUALITY 


1 1 tin t 



Unlike an untwisted elliptical planform wing, which has ellintir loadinn 

Educed less ?Ld?na e on?v al^V 'T' '° SS P ro P 8ller ha's minimi 
+^„ .1 i ! 1 9 °? ,Y , a+ 1+5 des 1 9 n advance ratio. An arbitrary propeller 
theory is needed to calculate its off-design point performance or the nerfor 

= + i ini tl v i ve s "h ,s ,? 

performanc e of a minimum induced los^ propel”; Shen appl IL?^' 9 " 
conditions and geometry calculated by the methods described before. 9 


-The axial and swirl components of the induced velocity at the blade ele- 

%n a"nu us n of Ihrll^V^ Cha " 9 ? S ° f a * ia ' a " d swfrl ^-ntum Ji+tVl 
given annu us of the slipstream equal to the axial and torque loadina of +h*> 

corresponding blade elements as shown * loading of the 


In figure 4: 


where (see figure 4) 


d_T 

dr 

= 2irrp V 

(u + a) 2FaV 



= P. y ^ 
2 

r u + aF 
[sin <j>j 

' Be ' 
2nr 4 

! 2irr C 

y 

(29) 

1 

r dr 

= 27rrpV 

(u + a) 2Ffira ' 



= P v 2 

'u + af 

' Be ) 

2nr C 

X 

(30) 


2 v 

sin 4>J 

2*r\ 


= tan - '( 

'_V_ iZj 
.ft r (1 - 

■ a) 

‘ a? > 

) 

(31 ) 

C ^ 

y 

cos 

4> - c d 

s i n 


(32) 

O 

X 

II 

c £ sin 

* + c d 

cos 

4> 

(33) 


In the absence of the propeller, the velocity in the flow field about th~ 
fuselage or nacelle is assumed to be given by an average axial component^ and 
an average radial component v at a distance r from the propeMer shaJ? 

We account for only the axial component 


u = u/V 


(34) 


—A— = 1 - C y 1 . a = Jc_ 

U + a 4 sin 2 <j> F 2irr 


(35) 


a C 


I - a* 


4 sintj) cos<)> F 


(36) 
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ii*/ * 


it* 


Equations simlliar to these appear In Glauert's article in Durand's "Aerodyna- 
mic Theory" (ref. 7), with the vortex spacing factor F in the numerator in- 
stead of the denominator, just as his widow and R. McKinnon Wood left them. 


The induced velocity components are evaluated at each radial station by an 
iterative process outlined below: 


At each value of £ the following are known: 


K, X, F, a, 3, u; c 0 = c. (a) , c, = c,(c„) 


Choose a 


d A 


Calculate <j) a = 0 - ct| 

Calculate c„, c. 

A d 


Calculate C , C (eqs. 32 and 33) 

y * 


Calculate a and a' (eqs. 35 and 36) 


Calculate <|> = tan 


-I 


I 


X ( u + a ) 


U (: - a') 


Calculate <J> - <f> 


a 


if * 0( - ^ > 0, a 2 < a. 


if 4> - <P- < 0, a 7 > a. 

a l a, 21 


Iterate until |<f> a - <f> a | is less than some small quantity. 


n 


Reta i n <b , C , C , 
r n' y ’ x ' 
'n n 


The wing theory analog of this computation is to suppose that the induced angh 
of attack at any spanwise station y of a non-e I I i pt ica I ly loaded wing of span . 
is given by 

, (c/b)c„ 


'induced 


/TT 


(37; 


(2y/b ) 

where c and c & are the cho rd and section lift coefficient at the same station. ! 

The quantity /] - (2y/b ) 2 vanishes at y = b/2 in the same way that F vanishes 
at r = R, and it may be shown that equation 37 yields «T nduced = C./ir(b 2 /S) for) 
an el I i pt ica I I y loaded wing of elliptic planform. i 


The values of <J>, C , C x , and a' are then integrated radially to find the 
thrust load and the power absorption of the propeller in the fuselage (or 
nacelle) flow field. These may be conveniently written in terms of coeffi- 
cients based on the shaft speed n (revolutions/sec) 


C T 


pn 2 D 4 


(D = 2R) 


(38/, 


292 


ORIGINAL PAGE IS 
OF POOR QUALITY 






ju 4, . 


k t r? hi 


r - r 

p 7TF 

pn D 






ll _ li fl - a '] 2 C 3 aC 


| d£ 4 [cos </) J y 

I H = III fj - a 1 ) 2 4 

j d£ 4 [cos (j) J ^ aC x (4 

t 

{^choice? 05 * Ch ° iCe ^ Preferred for numer ica I precision over the a, sin <fr. 


| The thrusting propeller is surrounded by a static pressure field wH-h =n 
jha<lestT mated* Its ZllV™' ^ downs+ream - Zoning (ref. 8) 


J?L = u- _= fi x/R 

P V 2 2 I 1 - , 

2 /(x/R) 2 + 


downstream; 
tractor prope I ler 


/(x/R)' 


upstream; 
pusher prope I ler 


Here x is distance downstream from the propeller. This axial pressure nred!cm+ 
causes the propel ler-beari ng-f use, age or nSdle to h/e a L^cy 


_ yX. Qr b 


buoyancy 5 ref 0 IfL/J [%J b 

and r s re,orence f ° r drag coefficients, 2 is the body length, 

”e net thrust^of y t he° S5 '“Vu® 3 at ths dist =" ce * b behind the body nose, 

me net thrust of the propeller-body combination is then given by. 


c T = c T 

n 


■^buoyancy 


while the installed efficiency of the propeller becomes 


(V/nD) 


.. Pi 9 ur ® 5 sh °ws an application of the arbitrary propeller theorv iust 8ps- 
cr.bed to the prediction of the performance of a scale SodeT ofTngh? Ii?! 
plane propeller when tested as an "isolated" propeller, and when runat the 

no th n a th7dT enta I/ e fUSe ' a9e - Th ' S ,heor V is computationally more Semand- 
ing than the design theory presented in the previous section since it reauirec 

.extensive estimates of the propeller airfoil section properties at several 
f radii, a good estimate of the three dimensional flow field surrounding the 


r i Ar/i lit - 


fuselage (or nacelle) at the propeller location, an Iteration procedure to de- 
termine the Induced velocities, and numerical integrations to determine blade 
loading and body buoyancy drag in the propeller pressure field. Limited ex- 
perience with it at M.l.T. shows that it gives reasonable results, and these 
are being experimentally confirmed (1979). In common with other radially 
graded momentum theories it fails to take account of the effect of circulation 
at every radial station on the downwash (or "inflow") at a particular station, 
but it is made to be consistent with the induced velocity pattern for a mini- 
mum induced loss propeller through Prandtl's analytic vortex spacing velocity 
fraction F rather than through tabulated values of Goldstein's circulation 
function. The next step up to a "prescribed" or "free" discrete vortex model 
of the "rotor" and its "wake" is much more diffcult. 


APPLICATIONS 


(I) Man powered airplane. Here we redesign the "Gossamer Condor" propeller. 
The design conditions, summarized in figure 6, correspond to climbing flight in 
ground effect at an angle of 1°; approximately 30 % of the 53.3 N (11.8 lbs) of 
thrust is required to overcome the component of airplane weight along the 
flight path. The figure shows the radial variation of profile D/L ratio and 
the radial gradients of the integrals l|, I 2 , J|, and J 2 . The design thrust 
coefficient, T c = 0.3175, requires a displacement velocity ratio, £ = 0.2671, 
which corresponds to a power coefficient, P c = 0.3914, and an efficiency, 
n = 0.8113. The powerplant output required is 328 watts (0.44 hp). 


Since the displacement velocity ratio is moderately large, it is worth- 
while to recalculate the thrust coefficient and the power coefficient using 
eqs. 21-24. The results are summarized in figure 7, which also compares the 
propeller geometry determined by the methods of this paper with the geometry 
actually employed. The agreement of blade angles is very good, especially 
when one takes into account the difference between the zero lift angles of the 
Clark Y airfoils assumed in the design calculations and the Stratford pressure 
recovery airfoils used on the "Gossamer Condor". In my opinion the propeller 
calculated here would be more efficient than the one actually flown. 


(2) Powered hang glider. Soarmaster, Inc. supplies a powerpack consisting of 
a West Bend (Chrysler) two stroke, single cylinder engine developing 7.46 kW 
(10 hp) at 10,000 rpm, a centrifugal clutch, a chain and sprocket reduction 
gear, and an extension shaft turning a pusher propeller. This is a suitable 
powerplant for hang gliders of 12 m (40 ft) span; figure 8 presents the options 
available for propellers intended to absorb the engine power at a flight speed 
of 13 m/sec (30 mph). The diameter of the direct drive propeller is limited 
to 690 mm (27 in) by a tip Mach number of 0.85; its efficiency is very poor 
because of the excessive disc loading. Gear reductions and larger propellers 
lead to progressive improvements in performance. Figure 9 gives the geometry 
of the largest propeller considered, a. 1372 mm (54 in) diameter propeller 
turned at 1946 rpm by a 9:37 sprocket pair driven at 8000 rpm. It has 617 mm 
(24 in) nominal pitch, and the typical wide root chord - narrow tip chord 
geometry of a propeller matched to a low advance ratio; this is in spite of a 
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design lift coefficient of I near the hub and 0.5 at the tip. Soarmaster 
supplies two propeller options: a 1067 x 483 mm (42 x 19 in) or an 1118 x 356 
mm (44 x 14 In) "laminar" airscrew, , both of fiber reinforced plastic. 

Table 2 summarizes the propeller parameters covered In this study: when 
two values are given for T c , P c , and r\, the second set corresponds to the 
■ improved velocity polygon geometry correspond i ng to eqs. 21-24; note the rela- 
J tlvely good agreement, even for £ values of more than I. Examination of 
I figure 9 and Table 2 suggests that still larger propellers and larger reduction 

( ratios would Improve climbing performance; this has to be balanced against the 
weight penalty and the reduction of ground clearance at the tail. 

^ (3) Motorsoarer . The Ryson ST- 1 00 is a 17.58 m (57.67 ft) span two seated 
| ai rcraft, with a flying mass of 748.4 kg (1650 lbs), fitted with a Hoffmann 
i HO-V-62 propeller of 1.7 m (67 in) diameter. This propeller has a low pitch 
I setting, a high pitch setting, and can be feathered for glider mode operation. 

1 Figure 10 shows three design points which might be considered in the selection 
i of such a propeller: sea level climbing performance at 40 m/sec (90 mph); 

? sea level top speed at about 68 m/sec (152 mph); and cruise at 15% power at 
Hull throttle at 1981 m (6500 ft) altitude and 65 m/sec (145 mph). The circled 
: points show the performance that may reasonably be expected from minimum 
^induced loss propellers designed for each of these flight conditions by the 
I methods of the paper. 

Figure II shows how a compromise propeller may be designed which will 
give nearly this performance at two of these points. The displacement veloci- 
ties are calculated assuming minimum induced loss loading and a somewhat pessi- 
mistic radial distribution of D/L ratio. Blade lift coefficients are assigned 
at £ = 0.3 and £ = 0.7 so that the blade chord to radius ratio, c/R, as given 
by equation 26, is the same for both flight conditions. The c/R ratio Is then 
calculated at other radii, assuming a linear radial variation of c^. 

Reasonable assumptions are then made about the radial variation of thickness 
to chord ratio, t/c, to give the radial variation of blade angle (eq. 28). 

The compromise c/R ratio and blade twist, A6, are then chosen to minimize 
differences between the two conditions. In general, highly loaded, low 
-advance ratio flight conditions demand high lift coefficients near the hub 
Betz (ref. 9) was of the opinion that Coriolis forces within the rotating 
blare boundary layer favored such a distribution. 
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TABLE I 

SYMBOLS AND NOTATION (follows Glauert; ref. 8) 

aV axial component of induced velocity (m/sec) 
a'ftr rotational (swirl) component of induced velocity (m/sec) 
B number of blades 

b wing span (m) 

c blade chord (m) 

cj section (profile) drag coefficient 

C£ section (profile) lift coefficient 

C[_ wing lift coefficient 

C p power coefficient (C p = P/pn-^D^) 

Cj thrust coefficient (Cj = T/pn^D^) 

C x blade element torque load coefficient 

Cy blade element thrust load coefficient 

D x drag; also propeller diameter (m) 

F slipstream velocity fraction (eq. 5) 

f vortex sheet spacing parameter (eq. 3) 

G circulation radial distribution function (eqs. 10,11) 

I 1 , 1 2 Lhrust loading integrals (eqs. 16,17) 

J | ,J2 power loading integrals (eqs. 18,19) 

L lift 

n revolutions per second 

P shaft power (kW) 
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power coefficient (P c = 2P/pV>jtr2) 

propeller shaft torque (Nm) 

propeller tip radius (m) 

propeller general radius (m) 

wing plan or fuselage cross section area (m 2 ) 

thrust (N) 

thrust coefficient (T c = 2 T/pV 2 ttR 2 ) 
axial velocity of fuselage flow field at r (m/sec) 
average axial velocity at r (m/sec) 
flight velocity (m/sec) 

radial velocity of fuselage flow field at r (m/sec) 

displacement velocity (m/sec); see fig. I 

resultant velocity at blade element^(W = V + s|r + w) (m/sec) 

induced velocity at blade element (w = a V + a ’fir) (m/sec ) 

slipstream velocity ( incrementa I ) (m/sec) 

velocity ratio (x = flr/V) 

spanwise location (m) 

section angle of attack (rad); a 0 (degrees) 

section blade setting angle (rad); 6° (degrees) 

circulation (m 2 /sec) 

displacement velocity ratio (£ = v'/V) 

efficiency (n = T c /P c = (V/nD)C T /C P ) 

advance ratio (A = V/QR) 

radius ratio (? = r/R) 

air density (kg/m^) 

blade solidity (o = Bc/2tit) 

helix angle (rad); cj> = $ - a 

shaft speed (rad/sec) 


TABLE 2 


POWERED HANG GLIDER PROPELLERS 
V = 13.41 m/sec (30 mph) 

7.457 kW (10 hp) @ 8000 engine rpm 
5 = 1.225 kg/m- 5 (760 mm Hg, I5°C) 


Gear 

Ratio 


9: 

:27 


: 27 

T\ 

37 


.690 


0 

2.5 
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F/ 6.4 PAD /ALLY GRADED AtOAiFNTUM THEORY 
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F/6 . 5 Radially graded momentum theory 

4 SCALE- MCCAUUY 1060 

- 1 


-105 

C p 


CAICULLATIOA/s H) 
fierce vombecr : 



0.1 A 02 


0 3 0.1 
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— V£T 
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FK*. 6 PROPELLER DES/GA/ FUJVCT/OUS 

MAR POU/£eSD A/RPLAVr 



0.021 

V/L 

o L 

o 



0.5 1 


R- 1. 905 *n ( 6.25 /?) 

Y * 5 m/sec (H2mM) 
ft* 1152 rad /sec (llO rptn) 

T - 53.3 V ( ll.e fhSj f°cli*<b) 
P * 1.170 Y 3 /tH* J 

L t * dli/dg Lf =1 2125 

/t = di z /df i 2 = o. osee 

J t = dj,/df J 4 = 13151 

J M *dJt/d£ J z * 0 5626 

T e - 0.3175 t-0 2611 

7> e * 0.3914<*. +! =■ 0.et13 

*^[/*|] * 0.6114 
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% t = 3.345— ~T~ l.fSf "I * 0.550 
J « 1.151 — ^ -3.146 *J * 0.578 
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